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Abstract 



In billiard systems with a flux line semiclassical approximations for the density of states 
contain contributions from periodic orbits as well as from diffractive orbits that are scattered on 
the flux line. We derive a semiclassical approximation for diffractive orbits that are scattered 
once on a flux line. This approximation is uniformly valid for all scattering angles. The 
diffractive contributions are necessary in order that semiclassical approximations are continuous 
if the position of the flux line is changed. 
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1 Introduction 



Billiard systems with flux lines are simple dynamical systems that can model the typical be- 
haviour of low dimensional quantum systems. One main application has been in the investiga- 
tion of universal properties of quantum systems that are chaotic in the classical limit 0. The 
high lying states of these systems have statistical properties that depend only on the symme- 
tries of the system and can be described by random matrix theory ^ ^ . Ordinary billiards 
are used to model systems that are invariant under time inversion and are described by the 
Gaussian orthogonal ensemble (GOE). The introduction of an additional flux line allows to 
study a further universality class since the flux line breaks the time-reversal symmetry. If the 
flux strength is sufficiently large (when considering levels in a fixed energy range) the corre- 
sponding universality class is that of the Gaussian unitary ensemble (GUE), but by varying the 
fiux strength one can investigate also the GOE-GUE transition or parametric correlations in 
the GUE regime. Another application is for integrable billiard systems where the introduction 
of a fiux line can break the integrability of the quantum billiard or to use billiards with a 
fiux line to model the properties of quantum dots (see e.g. 0). 

A major advantage in all these applications is that the classical trajectories are not changed 
by the presence of a fiux line, except for the set of measure zero that hit the fiux line. In a 
semiclassical analysis the same set of periodic orbits appears independent of the fiux strength 
and only the semiclassical contributions of these orbits are changed. For example, in approxi- 
mations to the density of states the periodic orbit contributions have an additional phase 27rma 
where a is related to the fiux strength and m is the number of times the orbit winds around the 
fiux line. However, this is not the only semiclassical effect of a flux line. Its presence leads also 
to wave diffraction that can be semiclassically described by an additional set of trajectories that 
are closed but not periodic and which start and end on the flux line. In semiclassical arguments 
these so-called diffractive orbits are often neglected since their semiclassical contribution is es- 
timated to be of order \fh, smaller than the contributions of periodic orbits. Such an estimate 
is, however, only valid if the scattering angle is not too close to the forward direction. In the 
forward direction diffractive orbits contribute in the same order as periodic orbits. 

In this article we investigate the semiclassical contributions of diffractive orbits to the density 
of states. We derive an approximation for orbits that are scattered once on a flux line. This 
approximation is valid for all scattering angles, also in the forward scattering direction. We 
show that these contributions are necessary in order to cancel discontinuities in periodic orbit 
contributions that occur when the position of the flux line is changed and crosses a periodic 
orbit. We also discuss the importance of diffractive orbits in the semiclassical limit. 

The article is organised as follows. In section 2 we consider the scattering on a flux line in a 
plane and derive a uniform approximation for the Green function. With this input we derive in 
section 3 a uniform approximation for semiclassical contributions to the density of states from 
isolated diffractive orbits that are scattered once on a flux line. As an example of a system 
with non-isolated diffractive orbits we consider in section 4 the integrable circular billiard with 
a flux fine in its center. The results are discussed in section 5. 

2 A flux line in a plane 

The scattering of a wave function on a flux line in a plane without boundaries has been studied 
in great detail since the problem was flrst treated by Aharonov and Bohm [0]. There exists, for 
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example, an exact integral representation for the propagator. A review on the quantum effects 
of electromagnetic fluxes is given in . In this section we review some of the results and derive 
a uniform approximation for the Green function that will be used in the following section. 

The Schrodinger equation for a particle with mass M and charge g in a magnetic field in 
two dimensions is given in Gaussian units by 
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m{r) = ^^(r) , (1) 
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and for a flux line the vector potential can be chosen as 
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VxA = ^6{r)z, V-A = 



(2) 



which describes a magnetic field with flux $ that is concentrated in form of a delta-function in 
the origin. 

In polar coordinates the Schrodinger equation with vector potential has the form 
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where a = q(^/{2'Khc) and k = \/2ME/h. The solutions that are regular in the origin are given 
by 
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where the normalisation is chosen such that 
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as follows from the orthogonality relations of the exponential and the Bessel functions 0. 
It is sufficient to consider only values < a < 0.5 since the solutions for other values of a 
can be obtained by a multiplication of a (^-dependent) phase factor and possibly a complex 
conjugation. 

With the normalised solutions (^) the propagator can be written down directly 
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and in the same way one obtains an exact representation for the Green function 
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where r< and r> are the smaller and larger values of r and tq, respectively. 

For further evaluations the representation in the basis of the angular momentum eigenstates 
is not convenient. One can derive, however, from the sum in (^) an exact integral representation 
for the propagator (see p]). It involves an integral over Hankel functions. In the following we 
use an approximation to this integral representation that is obtained after replacing the Hankel 
functions by their leading asymptotic form. This gives a semiclassical approximation for the 
propagator that is valid if < a < 1 and MrvQ/t ^ h, i. e. for not too long times or too close 
distances to the flux line. 
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The angles in (|^) have to be chosen such that |0 — 0o| ^ """j ^^nd the function K{z) is described 
below. The semiclassical propagator in (P) consists of two parts. The first part is almost 
identical to the propagator in a free plane with the difference of an additional phase proportional 
to a. Semiclassically it can be interpreted as the contribution from the direct path from Tq to 
r. The additional phase results from the dependence of the Lagrangian on the vector potential. 
The second term in (P) describes the scattering on the flux line and is discussed in more detail 
in connection with the Green function. Both terms are discontinuous in the forward direction 
10 ~ 00 1 = but the sum is continuous and uniformly valid for all values of and 0o. 
The function K{z) in (H) is a modified Fresnel function that is defined by 



K{z) = exp dy e'^" = -e"'"' erfc ie-'^'^z) , (9) 



and it has the following limiting properties 

^(0) = -; ^(^)~^' 1^1^^' -i<^^gW<T- 

An important alternative representation of K{z) is given by the integral z > 0) 

Due to this form the function K arises in semiclassical evaluations of oscillatory integrals in 
which a pole of the integrand is close to a stationary point. In figure ^ we show the real and 
imaginary parts of the function K{z) for positive z. The function has its largest absolute value 
at z = 0, and from approximately 2; = 3 on it agrees well with its asymptotic approximation 

With the relation between the propagator and the Green function in one can obtain a 
uniform approximation for the Green function. We consider the contributions from the two 
parts of the propagator in (Q) separately. These parts are called the geometrical and the 
diffractive part in the following 

G^ir, ro, E) ^ Ggir, r^, E) + Ga{r, r^, E) . (12) 



4 





Figure 1: Real and imaginary parts of the modified Fresnel function K{z) (full line) and its 
asymptotic approximation (dotted line). 



For the first part one can evaluate the integral in ([^ by stationary phase evaluation and one 
obtains, analogously to the propagator, the free semiclassical Green function modified by a 
phase 
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For the second part we express the Fresnel function by the integral (|Tl|) and arrive at 
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This double integral has a stationary point at {z,t) = {0,tci) where t^i = \/ Mj {2E){r + tq) is 
the classical time for the path from r to Vq via the origin at energy E. In order to evaluate the 
integral we expand the exponent in (|14D up to second order around the stationary point. Then 
the integral over t can be evaluated by stationary phase approximation, since the stationary 
point at t = tci is well separated from the pole at t = 0. In the integral over z, however, the 
z-dependence of the denominator has to be taken into account and this yields again a modified 
Fresnel function. In this way a uniform approximation for the diffractive part of the Green 
function is obtained 
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where again the angular coordinates have to be chosen such that |</) — (/)o| < vr. This ap- 
proximation for the Green function will be the ingredient for the derivation of semiclassical 
contributions of diffractive orbits in trace formulas in the next section. It is valid as long as kr, 
kro > 1. 

We consider now an approximation to this expression that is valid if |0 — 0o| is not too close 
to TT. Then the argument of the Fresnel function can be replaced by its leading asymptotic 



5 



term (|T0|) and yields 



Gd{r, ro, E) ~ ^ , , — 7^3x7- exp <^ ^fc(r + tq) + -(0 - 0o) + ^tt ^ (16) 



2Tikh^^JrrQ cos (^-^) I 2 

This approximation is of the general form that is obtained within the Geometrical Theory of 
Diffraction (GTD) (see e.g. fig) 



Gd{r, ro, ^) ~ ^ G,{r, 0, E) V{<f), 0o) Go(0, r-o, i?) . (17) 

In this theory the scattering is described by a free Green function from Vq to the scattering 
source at the origin, multiplied by a diffraction coefficient P(0, (po) that contains the information 
about the particular scattering process and a further free Green function from the origin to the 
point r. A comparison with (|16D shows that the diffraction coefficient for the present case is 
given by 

^l*^' '^0) = exp <^ — } . (18) 



cos 



(^) 



The term (|16|) can be interpreted as the contribution of a trajectory that runs from Vq to 
the origin and then to r. It is of order smaller than the contribution of the direct 

trajectory from ro to r in (|13|). However, the GTD-approximation breaks down in the forward 
direction |0 — 0o| = where the diffraction coefficient (|l^) diverges. This reflects the fact 
that the diffractive part of the Green function has a different leading asymptotic term in the 
forward direction. Here the diffractive trajectory contributes in the same order in k as the 
direct trajectory. The uniform approximation ([T5|) interpolates between these two different 
asymptotic regimes. 

In figure |^ we compare the different approximations to the Green function with the exact 
result (^ . The dotted hne is the geometrical part (0) . The approximation is already reasonably 
good, but one can still see a clear deviation from the exact curve. By adding the diffractive 
part in the GTD-approximation (dashed line) the difference becomes much smaller for most 
values of — 0o, except near the endpoints (j) — (po = ±tt where it diverges. This divergence 
is removed by the uniform approximation in (|15D (dashed-dotted line) for which the difference 
with the exact line can hardly be seen, even though we chose relatively small values of kr and 
kro- 



3 Semiclassical Contributions of Isolated Diffractive Or- 
bits 

For the derivation of semiclassical approximations in billiard systems it is convenient to apply 
the boundary integral method. It provides an alternative formulation of the quantum me- 
chanical eigenvalue problem in terms of an integral equation along the billiard boundary. All 
semiclassical contributions due to diffraction in billiard systems that go beyond GTD have so 



far been derived by this method |jTT 
iry 

(see e.g. p|, |1 
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The boundary integral method has been developed for billiard systems without internal 
fields 



but it can be extended to more general situations. In 15 Tiago et 
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Figure 2: Real and imaginary parts for different approximations of tlie Green function 
G{r,ro,E) with r = 1, tq = 2, = 5 and a = 0.4. Full line: exact result, dotted line: 
Gg{r,rQ, E), dashed line: GTD-approximation, and dashed-dotted line: uniform approxima- 
tion. 



al. derived an integral equation for billiard systems with magnetic field which is described by 
a vector potential in the Coulomb gauge. The resulting integral equation has the same form 
as for billiards without field. For a billiard system with a flux line and Dirichlet boundary 
conditions one obtains 

-2 [ ds'dnGo,{r',r,E)u{r') =u{r) (19) 

JdV 

where from here on we use units in which h = 2M = 1. The difference to the field free case 
is that the Green function in (|19]) is the one for a flux line in a plane and the solutions u{r) 
are the complex conjugate of the normal derivatives of the wave functions. The form of the 
integral equation in (|19]) is different from that in [^] since we applied an additional normal 



derivative to the equation in order to obtain a non-divergent integral kernel. Also the Green 
function differs by a factor (— Stt) from that in |jl5|. 

In ( P!^ the points r and r ' lie on the boundary, the position of r ' is parametrised by s\ and 
dfi denotes the outward normal derivative at r. The integral is evaluated along the boundary 
dV of the billiard system. The left-hand-side of (|TUp can be abbreviated by Qu{r) where Q 
is the integral operator acting on u. Equation (|19D is a Fredholm equation of the second kind 
and it has non-trivial solutions only if the determinant A(£^) := det(l — Q{E)) vanishes. This 
condition determines the quantum energies of the problem. For a summary of the properties of 
the Fredholm determinant A(i?) for two-dimensional billiards without fields with corresponding 



references see fiG 



From the condition of the vanishing of the Fredholm determinant one can obtain an expres- 
sion for the density of states in terms of the integral operator Q. It is given by 

Id 



d{k) = dUk) + -:77 Im V - Tr Q"(A;) , (20) 
71 dk ^-^ n 

n=l 



where 



Tr Q"(A;) = (-2)" / ds^ . . . dsn 5^,G«(r2, n, E) dn,G^{r^, rs, E) . . . dn^G^^r^ r„, E) . (21) 
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The sum in (^0[) is not convergent for real k, and for the derivation we assume that k has a 
sufficiently large imaginary part and that Re A; > 0. The density of states in (|20D has a smooth 
part and a sum over oscillatory integrals. A semiclassical evaluation of the oscillatory integrals 
yields the leading semiclassical contributions to the density of states from periodic orbits and 
diffractive orbits. These contributions can be separated by writing Ga as sum of its geometrical 
and its diffractive part. Then (|21|) consists of 2" terms. The periodic orbit contributions are 
contained in the term that contains only geometrical Green functions. By evaluating this term 
in stationary phase approximation one obtains the usual Gutzwiller expression for isolated 
periodic orbits Jl^, modified by an additional phase 27rma where m is the winding number of 
the orbit around the fiux line. The terms with I diffractive Green functions on the other hand 
contain the contributions of diffractive orbits that are scattered / times on the flux line. 

In this article we consider only diffractive orbits that are scattered once on a flux line. 
We thus have to evaluate an integral with a product of n — 1 geometrical Green functions 
and one diffractive Green function. This integral can be considerably simplifled by applying a 
composition law for the geometrical part of the Green functions. 



Gf{r, r',E)^ (-2)" / ds,... ds^ G,(ri, r E) d^^^Ggir^, r,, E) . . . dn^Gg{r, E) 

JdV 



(22) 



where the approximate sign denotes that the integrals are evaluated in stationary phase ap- 
proximation. The composition law was proved for ordinary billiard systems in [0 and it is a 
semiclassical version of the multiple reflection expansion of the Green function of Balian and 
Bloch. Since the presence of a flux line adds only a phase to the Green function which doesn't 
effect the stationary points in leading order, the composition law holds also in the present case 
and the phase is simply additive. The semiclassical expression for G^"^ is then (see W^ ) 



G^")(^2, r-i, E) = Y^ = exp <^ ikL - i-v - + ^a02i \ • (23) 

5„ A/STrfclMisI J 



Here 02i is the total winding angle of the trajectory around the flux line. It can be written as 
021 = 27rm + A0 where m is an integer and |A0| < tt. Furthermore, M12 is the 12-element of 
the stability matrix, L is the length of the trajectory, and i) is the number of conjugate points 
from Ti to r2- The stability matrix M is evaluated at unit energy and is energy independent. 
We use a tilde here in order to distinguish the quantities from the corresponding ones for the 
diffractive orbit. With (^21) one obtains the following expression for the partial contribution to 
the level density from diffractive orbits with one scattering event 

dpart = -lm^[ dsi ds2 dn,Gdir2, ri, E) dn^G'^^~^\ri, ra, E) . (24) 

TT afc Jgj, 

We insert now ( ^31) and (p!5|) with (0) for the semiclassical and the diffractive Green func- 
tion, respectively. The normal derivative gives in both cases in leading order a factor ik cos a 
where a is the angle between the normal and the outgoing trajectory. We obtain 

d A; cosai cosa2 sin(a7r) f .7r_ . , , 

dpart = Im — , ^ =- exp <^ -i-u + «a02i - > I , (25) 

27r3j|Mi2|(ri + r2) ^ ^ ^ ^ 
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where 



dsi ds2 
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The main part in the derivation consists in the evaluation of the diffraction integral /. If all 
three integrals are evaluated by a stationary phase approximation one obtains the contribution 
of the diffractive orbit in the GTD approximation. This expression diverges when | A0| = tt. In 
order to obtain a semiclassical approximation that is uniformly valid for all angles one has to 
take the dependence of the denominator on the integration variables into account. One can do 
this by expanding the exponent in (^Gf) up to second order in si, S2 and z, and the denominator 
up to first order in these variables. In the denominator one obtains in this way 



Acj) A0O , 1 
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This method corresponds to a particular choice of a uniform approximation. It yields an 
approximation that is correct at the 'optical boundary' |A0| = tt and in the GTD limit, and 
it interpolates between them. There are other possibilities to obtain a uniform approximation, 
for example by mapping the exponent onto a quadratic function in the three variables and 
then choosing an appropriate approximation for the amplitude function. This reflects the fact 
that uniform approximations are in general not unique. For example, a different interpolating 
approximation can be obtained by dropping the sine-term in (pTl) which corresponds to replacing 
it by its value at the optical boundary. 

The evaluation of the integral / is rather lengthy. It consists in carrying out the expansions 
and performing linear transformations of the variables such that, finally, the denominator de- 
pends only on one of the variables and the integral can be expressed in terms of the modified 
Fresnel function. A comparison with |]T^ shows, however, that the same diffraction integral 
occurs for the diffraction on billiard corners. For that reason, the calculations don't have to be 
done again and the result can be inferred from this paper. One has (c > 0) 
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(28) 

where b.c. denotes a boundary contribution whose origin is the discontinuity of the diffractive 
part of the Green function. It is cancelled by a corresponding boundary contribution from the 
geometrical part of the Green function. In (^) L and M denote the length and stability matrix 
of the orbit, respectively, and v is the number of conjugate points along the orbit, k is given 
by 



TrM-2 ^ ^ ' 
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TrM-2 ^ ^ ' 



(29) 
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and r = 1 — 2k. 

With ( ^81) and noting that the derivative gives in leading order a factor {iL) one obtains 
the final result for the contribution of a diffractive orbit ^ 
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(30) 



where fi = u + k, agrees with the usual definition of the Maslov index for periodic orbits. 

Formula ( PO] ) is the main result of this article. It describes the contribution of a diffractive 
orbit to the density of states and is valid for all scattering angles A0, and for < a < 1. A 
slightly simpler approximation can be obtained by dropping the sine-term in (^) and replacing 
the cotangent by a cosine, corresponding to the discussion after (pTf). 

The use of the full Fresnel function in (reffinal) is necessary in an angular region around the 
forward direction with width of order For angles outside this range the Fresnel function 

can be replaced by its leading asymptotic form which results in 
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(31) 



This agrees with the GTD approximation for the diffractive orbit [|18 



rf^(fc) = Re(-G'«(r2,ri)P(0i,02) 



(32) 



It is a contribution that is of order \/k smaller than the contribution of a isolated periodic 
orbit. In ( |52D G| denotes the contribution from the trajectory ^ to the geometrical part of the 
Green function where ri = r2 is the position of the flux line. 

In the opposite limiting case at \A(f)\ = n the uniform approximation is of the same order as 
that of a periodic orbit and it is discontinuous. As |0| goes through n the contribution makes 
a step of size 



dm 



A</>=7r-0 



dm 



A0=-7r+O 



2tL sin(a;7r) 
7r^|TrM-2 



sin (^kL — — yU + (2m + l)7ia^ . (33) 



In order to understand this discontinuity we consider an angle A(j) = tt—e which is inflnitesimally 
different from vr. By examining the linearised motion around the diffractive orbit one flnds that 
there is, in general, a periodic orbit in an inflnitesimal neighbourhood of the diffractive orbit. 
This orbit is obtained from the conditions that dq = dq' and dp = dp' — ke where the primed 
and the unprimed quantities are inflnitesimal perpendicular deviations from the starting point 
and endpoint of the diffractive orbit, respectively. From these conditions one flnds that the 
position of the periodic orbit is given by dq = —keMu/ [Tr M — 2). Depending on the sign 
of the left-hand-side, i.e. on the value of n, there can be two possible cases that are shown in 
flgure ^ for positive A0 near vr. 

If one moves the flux line in a way that A(f) goes through vr then the flux line crosses this 
periodic orbit at the same instant. For the periodic orbit the semiclassical contribution is also 
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Figure 3: Parts of a diffractive orbit and a nearby periodic orbit for k = and k = 1. 



discontinuous since the winding number of the periodic orbit changes. For a primitive orbit the 
winding number changes by +1 if k = and by —1 if k = 1. By writing the product of sines 
in (|3^) as a sum of two cosines one can show that the two discontinuities cancel exactly and 
the sum of both contributions is continuous. This shows that the uniform approximation for 
diffractive contributions is necessary in order to make semiclassical approximations continuous 
if the position of a flux line is changed. 



4 Diffractive orbits in the circular billiard 

A simple example in which diffractive orbits are not isolated is a circular billiard with a flux 
line in its centre. This system is integrable since the energy and the angular momentum around 
the centre are conserved. All diffractive orbits run from the flux line to the boundary and are 
reflected back directly onto the flux line. They appear in one- or more-parameter families, and 
their lengths are multiples of 2R where R is the radius of the circle. Contributions of these 



diffractive orbits to the density of states have been observed in |T9 |. 

We are interested in the leading order semiclassical contributions of the diffractive orbits. 
For this purpose one cannot apply formula (PD|) since the orbits have a stability matrix with 
trace two. Although one can in principle obtain the diffractive contributions from the boundary 
element method, it is now much more convenient to start from the torus-quantisation conditions 
for integrable systems (EBK-conditions). In this way one can obtain all leading order diffractive 
contributions to the density of states at once, even those for multiple diffraction. 

The exact solutions of the Schrodinger equation are given by the solutions of the flux line in 
a plane (§) with a different normalisation constant and the additional condition that the wave 
functions have to vanish at r = R. From this condition it follows that the energies are deter- 
mined by the zeros of the Bessel functions Em,n = ^^^mn/l^^), where km,n = j\m-a\,n+i/ R, 
and jj/^n is the n-th zero of the Bessel function with index u. 

The EBK-quantisation conditions yield semiclassical approximations to these energy levels. 
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They have the form 



dr Pr = 2Txh ( n + 



n = 0,1,2, 



(34) 



Due to the conservation of the angular momentum the first condition gives = hm and the 
second condition is evaluated with the energy conservation law E = {pl + {p^ — haY /r'^) / (2M) 
and yields 



\m — a\ arccos ■ 



\m — a\ 



T[ \ n + 



(35) 



which determines km,n as a function of the two quantum number m and n. 

From (|35| ) the periodic orbit contributions to the density of states are obtained by applying 
the Poisson summation formula 
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We consider now the main semiclassical contributions to the integrals in (|5B|). Altogether 
these are four contribution. One is the only non-oscillatory term M = N = 0, which yields the 
leading area-term for the mean density of states 



2 fkR 

dA{k) = — dm Vk^ 
TTfc Jo 



m^ 



IkR' 



(37) 



This agrees with the leading term in Weyl's law Ak/{2TT) with area A = nR^. 

The second main contributions arises from the stationary points of the integrals. A station- 
ary phase approximation gives the contributions of the periodic orbits that have been obtained 



also in 19 



[N/2] 



N=2 M+1 



Ak 



ttM 



irN 



N 



dpoik) = ^Y1 9m,n\I sin^ ^ cos (2NkRsm^ - + ^ ) cos(27rMa) (38) 



ttM 37r 



N 



where gM,N = 1 if M = A^/2 and gM,N = 2 if M 7^ A^/2. This factor arises from the fact that 
the stationary points for M = N/2 lie on the endpoint m = of the integrals in (|36| ) and give 
only half the contribution. 
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The remaining two contributions follow from the boundaries of the oscillatory integrals. 
They can be obtained by an integration by part which yields 

b. c. of I r dx g{x) e'^(") 1 = e'^^'^ , (39) 

for an upper end of an integration range. If z is on the lower end of an integration range the 
end point contribution is given by minus the right hand side of (pQ]). 
With f|39f) one obtains for the contribution from m = kr 



1 Vk^R"^ - m? exp {27iiM{m + a) - ^iN} 
dL{k) =2-^—7 /\/r o„-Ar„^„^^. m 



■nk 2-KiM - 2iN arccos ^ 

2R ^ sin (f AT) 

N=l 

_R 
~ 2" ■ 



\a — > —a\ 



m,—*kR 



(40) 



where the prime denotes that the term (A^, M) = (0, 0) is excluded from the sum since it 
corresponds to a non-oscillatory integral. The result in (|40|) can be identified with the perimeter 
term in the asymptotic expansion of the mean density of states. 
From the other end points of the integrals at m = one obtains 



dd{k) = TTT T7 N \ ^^^iMa + 2iNkR - -^iN \ + [a 



M,N=-oo 2 

= dlik) + d^k) . (41) 

Here the terms M = N/2 have to be excluded from the sum since they have already been 
taken into account by the stationary phase evaluation. For convenience the sum is split into 
two parts, corresponding to a summation over even and odd values of N, respectively. 

For the first part is replaced by 2A^, the summation index M is shifted by N, and terms 
for positive and negative values of M are combined 



\a —a\ 



R sin[27raMl f ,w , ^ Svr, 

N=-oo M=l ^ 

= ^(1 - 2a) ^ sin (2N{2kR - sm{7ia2N) , (42) 

^ N=l ^ ^ 

and for the second part is replaced by 2A^ + 1, the summation index M is shifted by A^, and 
terms for M > 1 and M < 1 are combined 

2R ^ ^ sinf7ra(2M - 1)1 f , ,^ ,^ SttJ , 

d2k) = -^ J2 E 2M-1 '^ expU2N + l){na + 2kR-—)[ + [a^-a] 

N=-oD M=l ^ ^ 

= — sin ( {2N + l){2kR - y ) j sin(7ra(2Ar + 1)) . (43) 
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The formulas for evaluating the summations in ( pO]) 
valid for < a < 1. 

Collecting all results, the complete trace formula is given by 



and (^31) can be found in |]2T[ and are 



oo [A^/2] 



d{k) « d{k) + 2^ 2^ 

N=2 M=l 



Ak , ^ ttM / , ttM 
W siir — — cos 2NkR sin ■ 



+ ^ ) cos(27rMa) 



vriV iV V N 

2R / Stt \ °° 2R 

+ V — sin 2A^fci? N sin(7rA^a) - 2a V — sin {4NkR - SttN) sm{27rNa) . 

^—^ 71 \ 2 1 ^ — ^ TT 

N=l ^ ^ N=l 

(44) 

It consists of the mean level density, the contributions of periodic orbits, and the contributions 
of diffractive orbits which are by an order smaller than those of the periodic orbits. The 



formula is similar to the corresponding one for a harmonic oscillator with a flux line p2 
The diffractive orbit term in 



can be given a more direct interpretation since it can 
be transformed into a sum of delta functions by using the Poisson summation formula. One 
obtains 



ddik) 



R oo 
It ^ — \ 

2^ ^ 



^ , kR 3 a 



n=— OO 
oo 



n 



^ , kR 3 a 

-6 { hn 

n 4 2 



aR \ ^ 

TT ^ 



l2kR 



- + a + n 



2kR 



TT 



-- a + n 



(45) 



The argument of one of the delta functions can be identified with the semiclassical quantisation 
condition for eigenvalues with vanishing angular momentum (m = 0): kR ~ 7ra/2 + 37r/4 + ?7,7r. 
Since the full trace formula produces delta peaks at the eigenvalues given by the EBK conditions 



(|35|), the diffractive orbits have the following role: they contribute to peaks at eigenvalues with 
vanishing angular momentum and they cancel wrong peaks that are produced by the periodic 
orbit sum. 

The fact that the periodic orbits produce wrong peaks can be understood by noting that 
the periodic orbit contributions (^) are invariant under a —a. Due to this symmetry the 
periodic orbits give also rise to wrong peaks at kR ?a — 7ra/2 + 37r/4 + rnr. The first two terms 
in (^) provide a correction to this failure in the approximation for states with zero angular 
momentum. The other two terms in (^Sf ) are necessary in order that the approximation is 
invariant under a ^ 1 — a, since the spectrum is invariant under this replacement. 

A test of the trace formula (|4^) is presented in figure ^ which shows a comparison between 
quantum mechanical and semiclassical results. The value of a is chosen to be a = 0.25, since 
for this value the periodic orbit contributions start at L = 8R. The difference between the two 
curves cannot be seen on this scale. 



5 Conclusions 

One general property of uniform approximations for isolated diffractive orbits is that they make 
semiclassical approximations continuous as a system parameter is changed. They cancel discon- 
tinuities in semiclassical contributions of periodic orbits. In billiards with concave boundaries 
or with corners those discontinuities are connected with the appearance or disappearance of 
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Figure 4: Fourier transform of the oscillatory part of the density of states for the circular 
billiard with a = 0.25 and R = 1 (full line) in comparison with the semiclassical approximation 
(dotted line). Inset: Magnification of one peak. 

periodic orbits. In the present case the discontinuity is due to a phase change in periodic orbit 
contributions. Besides that, the diffraction on a flux line is very similar to the diffraction on 
corners. The final formula is expressed in terms of the same interpolating function, but it is 
also simpler since it consists only of one term instead of four. This suggests, for example, that a 
semiclassical study of two-dimensional systems which are non-intcgrabic but in which the clas- 
sical motion is still restricted to a two-dimensional surface in phase space might be performed 
more easily on billiards with a flux line instead of pseudo-integrable polygonal billiards. 

In this article we considered only diffr active orbits that are scattered once on a flux line. The 
same method can be used to obtain semiclassical contributions of diffractive orbits with multiple 
scattering, but the formulas become increasingly more complex and have to be expressed in 
terms of multiple Fresnel integrals. The treatment of an integrable billiard with a flux line like 
the circular biUiard is much simpler. There all leading order diffractive contributions including 
those from multiple diffraction can be obtained in one step from the EBK conditions. Moreover, 
they can be summed up and shown to contribute only to states with zero angular momentum. 

One remaining question is whether in semiclassical arguments concerning spectral statistics 
in chaotic systems it is sufficient to consider only periodic orbits. One argument for this is that 
the semiclassical contributions of most diffractive orbits (GTD region) are by an order 1 / \/k 
smaller than those of periodic orbits where k is the wavenumber. For scattering angles near 
the forward direction they contribute more strongly, up to the order of a periodic orbit, but 
the corresponding angular regime of this transitional region decreases proportional to l/\/k. 
So one might argue that diffractive orbits become less and less important in the semiclassical 
regime. 

Let us discuss this point in more detail. By applying the trace formula in order to resolve 
adjacent energy levels with wave numbers of order k one has to take into account all orbits 
up to the Heisenberg length Lh (x k. In a chaotic system the number of diffractive orbits 
increases exponentially with the orbit length. If we assume that the scattering angle is uniformly 
distributed then the relative number of diffractive orbits in the transitional region decreases 
hke 1/Vk. However, the total number of diffractive orbits in the transitional region is still 
increasing exponentially due to the exponential increase of all diffractive orbits. A similar 
argument can be applied to orbits with multiple scattering. It is not obvious that these orbits 
can be neglected. Moreover one can show that even in the GTD approximation diffractive 
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orbits have a non- vanishing influence on spectral statistics in the semiclassical limit 
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